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ABSTRACT 

A system  of  approximate  equations  for  the  determination  of  thermal 
stresses  in  piezoelectric  plates  with  large  thin  films  of  a different 
material  plated  on  the  surfaces  is  derived.  The  plate  equations  jure  ob- 
tained by  making  a suitable  expansion  of  the  pertinent  variables  in  the 
thickness  coordinate,  inserting  the  expansion  in  the  appropriate  varia- 
tional principle  and  integrating  with  respect  to  the  thickness  in  the 
manner  of  Mindlin.  Conditions  resulting  in  both  extensional  and  flexural 
stresses  are  considered  and  the  full  anisotropy  of  the  quartz  is  included 
in  the  treatment.  The  particular  case  of  purely  extensional  thermal  stresses 
resulting  from  large  electrodes  of  equal  thickness  plated  on  the  major 
surfaces  of  doubly-rotated  quartz  thickness-mode  resonators  is  treated  in 
detail.  The  changes  in  resonant  frequency  resulting  from  the  thermally 
induced  biasing  stresses  and  strains  are  determined  from  an  existing 
perturbation  equation.  Calculations , using  the  newly  defined  first  temp- 
erature derivatives  of  the  fundamental  elastic  constants  of  quartz, are 
performed  for  large  gold  electrodes  on  doubly-rotated 


quartz  plates. 

) 


1.  Introduction 

A perturbation  analysis  of  the  linear  electroelastic  equations  for  small 
fields  superposed  on  a bias  has  been  performed*.  The  change  in  resonant  fre- 
quency due  to  any  bias  such  as,  e.g.,  a residual  stress  may  readily  be  ob- 
tained from  the  resulting  equation  for  the  first  perturbation  of  the  eigen- 
value if  the  bias  is  known,  the  use  of  this  perturbation  equation  has  already 

been  shown  to  be  extremely  accurate  in  the  determination  of  changes  in  the 

2 

surface  wave  velocity  of  crystals  due  to  flexural  biasing  stresses  . 

In  this  paper  a system  of  approximate  plate  equations  for  the  determina- 
tion of  thermal  stresses  in  thin  piezoelectric  plates  coated  with  much  thin- 
ner films  is  derived  in  the  manner  of  Mindlin3  5.  The  resulting  approximate 
equations  simplify  the  treatment  of  many  thermal  stress  problems  considerably, 
and  the  three-dimensional  detail  not  included  in  the  approximate  description 
is  not  deemed  to  be  important  for  our  purposes.  In  order  to  keep  the  deriva- 
tion of  the  thermoelastic  plate  equations  clear  and  not  introduce  extraneous 
complications  that  can  lead  to  confusion  at  the  outset,  we  first  ignore  the 
elastic  constants  that  cause  coupling  between  shear  and  extension  in  the 
constitutive  equations.  Both  extensional  and  flexural  plate  equations  are 
obtained.  After  we  obtain  the  thermoelastic  plate  equations  under  the  afore- 
mentioned simplifying  assumptions,  we  extend  them  to  the  general  anisotropic 
case.  It  should  be  noted  that  coupling  between  shear  and  extension  exists  and 
is  important  even  in  the  case  of  rotated  Y-cut  quartz.  The  general  aniso- 
tropic version  of  the  thermoelastic  plate  equations  is  applied  to  the  case  of 
purely  extensional  thermal  stresses  arising  from  large  rectangular  identical 
electrodes  on  doubly-rotated  quartz  plates  and  a simple  exact  solution  of  the 
extensional  plate  equations  is  obtained.  This  solution  of  the  static 


approximate  thermoelastic  plate  equations  for  the  arbitrarily  anisotropic 
plated  crystal  plate  with  rectangular  electrodes  is  readily  shown  to  hold 
for  large  electrodes  of  essentially  arbitrary  shape.  Hie  approximate  three- 
dimensional  displacement  field  resulting  from  the  solution  of  the  plate 
equations  is  readily  determined  from  the  description.  This  three-dimensional 
displacement  field  is  required  in  order  to  obtain  the  change  in  resonant 
frequency  due  to  the  thermally  induced  biasing  state  from  the  equation  for 
the  first  perturbation  of  the  eigenvalue. 

Since  the  total  change  in  frequency  of  quartz  plates  is  due  not  only  to 
the  thermal  stresses  and  strains  but  also  to  the  change  in  the  fundamental 
elastic  constants  with  temperature,  in  order  to  complete  the  above-mentioned 
calculations  of  the  resulting  actual  change  in  frequency  per  degree  change 
in  temperature,  the  newly  defined**  temperature  derivatives  of  the  fundamental 
elastic  constants  of  quartz  must  be  employed.  These  newly  defined**  temper- 
ature derivatives  of  the  fundamental  elastic  constants  of  quartz  are  used, 
along  with  the  aforementioned  thermal  stress  analysis  of  the  quartz  plate 
and  the  perturbation  equation,  in  performing  calculations  of  the  change  in 
frequency  per  degree  change  in  temperature  for  rectangular  gold  electrodes 
on  doubly  rotated  quartz  plates.  TSie  results  indicate  that  there  are  whole 
ranges  of  doubly-rotated  orientations  for  which  the  change  in  frequency  with 
temperature  due  to  the  electrodes  is  more  than  an  order  of  magnitude  less 


than  that  of  the  AT-cut 


T 


3. 


2.  Perturbation  Equations 

Die  equation  for  the  first  perturbation  of  the  eigenvalue  obtained 
from  the  perturbation  analysis  mentioned  in  the  Introduction  may  be 
written  in  the  form 

\ m v2v  ® = t ■ v (2,l) 

where  and  cu  are  the  unperturned  and  perturbed  eigenfrequencies, 
respectively,  and 

■V  ■ I ds  + J dv,  (2.2) 

in  which  S is  a surface  enclosing  a volume  V.  In  (2.2)  N_  denotes  the 

L 

unit  normal  to  the  undeformed  surface  at  the  reference  temperature,  u^ 
denotes  the  mechanical  displacement  vector  and  cp  denotes  the  electric 

~JL  ~jj 

potential.  Die  linear  stress  and  electric  displacement  vector 
are  given  by  the  usual  linear  piezoelectric  constitutive  relations 


V " 2bYMQfU Of,  M + eMLYCf>,M  ’ 


3 « e u - e cp 
L LMy  Y,M  LMy,M  * 


(2.3) 


and  as  usual  satisfy 


V, 


Pu, 


„ If 


L Y L,L 


0. 


(2.4) 


Die  quantities  c e and  e denote  the  second  order  elastic, 

2 LyMCr  MLy  LM 

piezoelectric  and  dielectric  constants,  respectively,  and  p denotes  the 
mass  density.  Die  nonlinear  contributions  to  the  Piola-Kirchhoff  stress 
tensor  and  reference  electric  displacement  vector  take  the  respective 
forms1' 
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\y  “ ^Ly mu  + ^lymq')uo',  m + ®mly^,  m ’ 

Jr  m e u - € cd  (2  5) 

L LMY  V,M  LM  , M , 1 * 1 

Where  c ^ and  are  effective  constants  that  depend  on  the  biasing 

state'*'  and  Ac  denotes  a small  change  in  the  fundamental  elastic  constants 

due  in  this  instance  to  a change  in  temperature.  Clearly,  the  total  dynamic 
Piola-Kirchhoff  stress  tensor  and  reference  electric  displacement  vector 


are  given  by 


K - + k"  3 - I£  + 5^  . 

LY  LY  lX  * L L L 


The  vector  g^  denotes  the  normalized  mechanical  displacement  for  the  nth 
unperturbed  mode  and  5^  denotes  the  normalized  electric  potential  for  the 


lith  mode, 


where 


ur  -41 

g »*  . Jt  „ 2L 

gY  N » N * 

n n 

Nn  " J pu?^  dv* 


n J y Y 

The  normalized  stress  tensor  k*^  and  electric  displacement  vector  d^,  both 

uy  L7 

for  the  normalized  nth  mode,  are  given  by 

# - c a*  + e & 

LY  LYMCr  Of,  M MLY  , M * 

^ " *LMYgY,M  ” ®LM^M  * (2,9) 

The  upper  cycle  notation  for  many  dynamic  variables  and  the  capital  Latin  and 
lower  case  Greek  index  notation  is  being  employed  for  consistency  with  Ref. 1, 
as  is  the  remainder  of  the  notation  in  this  section.  Hie  fact  that  the  capital 
Latin  and  lower  case  Greek  indices  refer  to  the  reference  and  intermediate 
positions  of  material  points  respectively,  is  not  important  here,  and  in  this 
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work  they  may  be  used  interchangeably.  We  employ  Cartesian  tensor  nota- 
tion, the  summation  convention  for  repeated  tensor  indices,  the  convention 
that  a comma  followed  by  an  index  denotes  partial  differentiation  with 
respect  to  a reference  coordinate  and  the  dot  notation  for  differentiation 
with  respect  to  time. 

When  the  electrical  and  electroelastic  nonlinearities  are  ignored  we 

have 


eLM  = ° * ®MLY  ° ’ 


(2.10) 


and 


„ 1 , 1 

c = T 6 + c E.  + c w +c  w . 

LyMo  LM  yor  3LYM0AB  AB  2l*YKM  Of,  K Y, K 


(2.11) 


where 


TLM-  2lmkNEkN  ' ^ 2 (WK,N  + WN,K)  * 


(2.12) 


c,  denote  the  third  order  elastic  constants  and  w denotes  the  static 

jLyMOrAB  K 


biasing  displacement  field.  Unis,  in  this  description  the  present  posi- 
tion £ is  related  to  the  reference  position  X by 

y(X  , t)  = X + w(X.)  + u (X  , t)  . (2.13) 

£ Ii  ~ ~ L ~ L 


For  the  electroded  crystal  plate  with  shorted  electrodes.  Which  is 

7 


of  interest  here,  the  boundary  conditions  take  the  form 

NAY  - Vh/Y^BCDUC,DA-2h,p/iV'  S’*0'  °n  S»  (2*14) 

where  A,  B,  C,  D take  values  confined  to  the  surface  of  the  crystal  plate 
and  skip  the  value  associated  with  the  normal,  2h ' is  the  thickness  of  the 
electrode  and  p ' the  electrode  mass  density.  Since  the  electrodes  are 
isotropic,  the  plating  stiffnesses  are  given  by 

(2.15) 


Y»__,  ” X,6.  6 +ji ' (6__6  _ + 6__6  ) , 

ABCD  O AB  CD  AC  BD  AD  BC  ’ 


Where  the  plate  Lame  constant  X/  is  given  by 


f 


r 


6. 


X7  - + 2\i')  , (2.16) 

o 

and  X7  and  p7  are  the  Lame  constants  for  the  plating  material,  the  boundary 
conditions  satisfied  by  the  normalised  pth  eigensolution  of  the  unperturbed 
linear  system  are  given  by 

VV  Vh"'^CD4,D-2h''1'^-  »"S-  l2-17> 

Taking  u^  and  cp  to  be  the  unperturbed  normalized  pth  eigensolution  g^ 
and  f^,  respectively,  substituting  from  (2.6)p  (2.14),  (2.17),  (2.5)2 
and  (2.10)  into  (2.2)  and  employing  the  divergence  theorem,  we  obtain 


H = - [ k"  g1*  dV . 

U J T\/  T. 


J Ly  y,  i* 
v 


(2.18) 


Now,  the  substitution  of  (2.5)^  and  (2.10)2  along  with  the  aforementioned 
selection  of  u and  $ into  (2.18)  yields 


” J ^LyMOf + A2LyMQf^9Qf,  M^y,  l dV  ’ 


(2.19) 


Which,  with  (2.1),  (2.11)  and  the  given  small  change  in  the  fundamental 

elastic  constants  Ac.  ..  due  in  this  instance  to  a change  in  temperature^, 
2LyMof  ’ 

gives  the  change  in  resonant  frequency  due  to  a biasing  deformation  and 
the  change  in  elastic  constants  resulting  from  a change  in  temperature. 


3.  Static  Plate  Equations 

A schematic  diagram  of  the  plated  crystal  plate  is  shown  in  Fig.l 

along  with  the  associated  coordinate  system.  Referred  to  this  coordinate 

3 

system  the  static  form  of  Mindlin’ s plate  equations  may  be  written 


K(n)  n1f(n'1)+P(n)  = n _ _ n i , 
KM,A  *2J  PJ  °»  n °>1>2  > 

where  A,  B,  C,  D take  the  values  1 and  3 and  skip  2 and 
K£>  ' 1 ❖iJ  aX2  > Pj”’  ' • 


The  linear  thermoelastic  constitutive  equations  for  the  mth  order  stress 


resultants  for  homogeneous  temperature  excursions  (T  - T ) take  the  form 

o 

2 

K*®)  = c y H E^")  - d(m)vT_(T-T  ) , (3.3 

IJ  IJKL  i—j  mn  KL  IJ  o 

n»0 


IJ  o 


where  we  have  taken  the  liberty  of  dropping  the  lower  script  2 on  the 

c . the  v„  denote  the  thermoelastic  constants  and 

IJKL’  IJ 

H = 2h  (m+n+1 (m  + n + 1)  m + n even 
mn 

0 m + n odd  , (3 

d ^ = 2h^n+D/(m  + 1)  m even 

0 m odd  . (3 


The  nth  order  plate  strains  take  the  usual  form 

C "I  <^<^^  + 1>^62LW^1)+62KWLn+1,>Ip 

where  2 3 

w - Y Xnw(n)  w - Y Xnw(n) 

A L 2 A » 2 Z 2 2 » 


2 3 

E - V vnE(n)  E « V vnE(n) 

BAB  L 2EAB  » E22  L 2^22  * 

n-0  n*0 


(3.8) 


The  thermoelastic  constants  v are  related  to  the  coefficients  of  linear 

IJ 

g 

expansion  or  by  the  usual  relation 

Kli 


:ijkl“kl 


f 
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The  static  form  of  Mindlin1  s simplified  extensional  equations  for  the 
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very  thin  electrode  platings  may  be  written 

Or'-°>  pbo)/-^b]!v  > 


where 


-2hV  - 

AB  ABCD  CD 


E-  -2hV^(T-To), 


(3.10) 


(3.11) 


and  since  the  platings  are  isotropic  in  the  plane 


Y ' = \‘  6 6 Lt ' (6  6 +6  8 


'o''AB”CD  ■"  ' AC  BD  AD  BC)  , 

+ 2n')  , V*'  = [2tl'(3X/+2tl/)/(X'  + 2lA,)]a'6 


AB 


(3.12) 


Which  are  for  the  upper  plating  and  similar  equations  hold  for  the  bottom 
plating,  but  with  double  primes  replacing  the  primes.  In  order  to  obtain 
the  static  equations  for  the  plated  crystal  plate  we  need  the  boundary 
conditions  at  all  interfaces.  The  traction  conditions  take  the  form 

K'o/J-O,  k!  (-h#)-0, 


2A 

H //  , 


K2A(“h  ) = K2A(h)  ’ ^A01  } = K2A("h)  ’ (3.13) 

and  the  conditions  of  continuity  of  mechanical  displacement  take  the  form 

.«»'  (0)  .v..(D  ,.2  (2) 


w. 


WA  +hWA  +hWA  > 


(0)  (0)  . (1)..2  (2) 
w = w - hw  + h w„  , 
A A A A ’ 


(3.14) 


.(0)'  (0)' 


since  the  thickness  displacements  of  the  thin  platings  w^  and  w^  ' take 
place  freely. 

In  this  section  in  order  to  keep  the  derivation  of  the  extensional 
and  flexural  equations  for  the  plated  crystal  plate  clear  we  simplify  the 
treatment  somewhat  and  ignore  the  elastic  and  thermoelastic  constants  that 
cause  coupling  between  shear  and  extension  in  the  constitutive  equations. 
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This  sinplification  and  attendant  restriction  of  the  equations  is  removed 
in  the  next  section.  When  flexure  and  extension  are  uncoupled,  from  (3.1) 
we  may  write  the  extensional  plate  equations  in  the  form 

1S£,)A-K22)+P21)  = °> 


■jEJa-  2K2B)+FB2)-°* 


(3.15) 


In  order  to  allow  for  the  free10  thickness  strains  and  E.!2^ , we  take 

22  22  * 


c(0) 

22 


.(2) 


0 , K'~  - 0 . 

> 22 


(3.16) 


In  addition,  in  order  to  eliminate  the  first  order  extensional  equation 
in  (3.15)  completely11,  we  take 

(3.17) 

Prom  (3.16)  and  the  reduced  form  of  (3.3),  with  (3.4)  and  (3.5),  we  obtain 


_.(0)  C21  _(0)  C23  ,,(0)  V22  m , 

.<2).  °21  -(2)  C23  (2) 

22  °22  11  33  ’ 

the  substitution  of  Which  in  the  reduced  form  of  (3.3)  yields 
.*  .(0)  . * .(0)i  . 1 .3  * _(2)  * _ (2) 


(3.18) 


,(0) 

'll 


K„.  *2h[c11E11  +c13E33  ] +3h  [Cj.jEj.j_  + c13E33  1 - 2hVn(T-  Tq)  , 


*iV  “ 2htCl3EU)  +C33E33>]  + 1 +C33E33)]  ’ ^33  (T  ‘ To)  > 

■S’  -!  h3[c*iAV  ***iAV' + ! h5[cuEu) +c33b33)]  -fbMi^-v » 

-!  h3[c*!3EU)  + ! h5[Ci3EU>  +C33EE)]  *!  ^33^- V » 

K(o)  - 2ho  rp<o)  + id  e(2)i 

K13  2hc55LE13  3 ®13  J * 

K(2)  - 2h3e  FI  P{0)  + b i -(2)"|  ,, 

K13  2h  C5s|_3  *13  + 5 ®13  J » (3.19) 


rr 


10. 


where 


CU-°11 


'12 


'22 


'13  13 


c c 

-mi  . C*  -c 


'32 


'22 


33  33  c 


22 


VU“V11-V22C12/C22>  V* 


33  V33  " V22C32/C22  * 


(3.20) 


From  (3. 15), with  (3.16)  and  (3.17),  we  obtain  the  resulting  second  order 
extensional  equations 


.(0)  . „<0) 


(2)  . „ (2) 


V + PB  -°'  KAB,A  + 


K =o. 


B 


(3.21) 


, (n) 


In  this  approximation  from  (3.6)  the  pertinent  plate  strains  EftB  (n=0, 2) 
take  the  form 

^ tn\  ! n \ Its  \ Its  \ 1 Its  \ lrs\ 

(3.22) 


E(n)=w(n)  E(n)-w(n)  E(n)  = i (w(n) +w(n)) 

E11  *1,1 » E33  _W3,3»  E13  2 (W1,3+W3,1)  * 


When  written  out, the  constitutive  equations,  (3. 11), for  the  upper  electrode 


plating  take  the  form 

,(0)  ' / w. <o> ' , /.. (0) _ /..*/ 


07'  -2h,[(X'+2li,)w''';  +l*w'  ' ] - 2h  v (T-  T ) , 
11  O 1,1  O 3,3  O' 


K^)/-2h,[(X'+2ji,)w(°^  +XV(0)  ] -2h#V*'(T-T  ) , 
33  O 3,3  O 1,1  O 


4V'  - »v«*£}/+,£>'). 


(3.23) 


and  similar  equations  exist  for  the  lower  electrode  plating,  but  with 
primes  replaced  by  double  primes.  In  the  case  of  these  extensional  equa- 
tions  the  displacements  w^0^  and  w^°^  , in  the  upper  and  lower  electrode 
platings,  respectively,  are  given  by 

.(0)'  (0)  . . 2 (2)  (0)'_  (0)  v2„(2) 


wJ  ’ -w'-'  +h-w'“'  , w-'  -w'“'+hw'  . 

A A A ' A A A 


(3.24) 

The  equations  for  the  plated  crystal  plate  are  obtained  by  employing  (3.13) 
in  (3.10)  for  both  the  upper  and  lower  platings  and  then  inserting  (3.10) 
for  both  platings  in  (3.2)2,  with  the  result 


1 


11. 


+ *£,a.  fb2)-i'2<,)a  +*S,,a>- 


(3.25) 


The  substitution  of  (3.25)  into  (3.21)  yields 

[K^)+KM)/+C\a-°>  I>fi)+h2(^2)/+^)l  ,A-°»  {3-26) 

Which  are  the  second  order  extensional  equations  of  equilibrium  of  the 
plated  crystal  plate.  Clearly,  Eqs. (3.26)  are  consistent  with  the  inte- 


gral forms 


where 


f v Ca,-°>  kCa-°> 


(3.27) 


V(0) -K(0)  + K(0)'+K(0)"  V<2)_K<2)  2,(0)'  (0)' 

*AB  KAB  +KAB  +KAB  » *AB  KAB  +h  (KAB  + KAB  ]>  (3'28) 

and  Nft  denotes  the  outwardly  directed  unit  normal  to  the  undeformed 
position  of  the  closed  curve  c in  the  plane  of  the  plate.  Hence,  on  an 
edge  of  the  plated  crystal  plate,  the  traction  boundary  conditions  that 


accompany  (3.26)  are 


rf  - c 

»aC  - k2) 


(3.29) 


where  the  represent  applied  traction  resultants  on  the  edge.  The 

alternative  displacement  conditions  are  on  w^°^  and  w^2\  The  specif ica- 
tion  of  appropriate  combinations  of  these  prescribed  conditions  may 
readily  be  shown  to  be  unique  to  within  homogeneous  rigid  plate  rotations 
of  zero  and  second  order  by  means  of  the  usual  Neumann  type  procedure12’ 13,4t 
At  an  edge  of  discontinuity  separating  one  region  from  another,  we  require 


the  continuity  of 


N jv(0)  „v(2)  w(0)  w(2) 

NA*AB  » "a*AB  • B » B * 


(3.30) 


12 


Since  the  upper  and  lower  electrodes  can  have  different  thicknesses 
2h ' and  2h*,  flexure  can  occur  even  though  the  temperature  change  (T  - Tq) 
is  homogeneous.  In  order  to  obtain  the  flexural  equations  we  first  write 
the  plate  equations  in  the  form 


Km!a  + P20)-0>  Km!a-^B,+PB1)  = °*  K-,.-2K,(i)+pi2)-0. 


* A2,A 


22 


(3.31) 


In  order  to  allow  for  the  free  thickness-strains10,  we  take 


_U) 

22 


0 . 


(3.32) 


In  order  to  eliminate  the  second  order  flexural  equation  in  (3.31)  com- 
14 

pletely  , we  take 

K^-O.  (3.33) 


In  addition,  as  usual  in  the  elementary  theory  of  flexure  we  take  the 
thickness-shear  plate  strains  E^0*  to  vanish15.  Which  with  (3.6)  yields 


,(°> 

2, A * 


(3.34) 


Prom  (3.32)  and  the  reduced  form  of  (3.3),  with  (3.4)  and  (3.5),  we  obtain 

c. 


E. 


d)  lil  (1)  ~23  _(1) 

22  -c22EU  -c22E33  * 


(3.35) 


the  substitution  of  Which  in  the  reduced  form  of  (3.3)  yields  the  flexural 
constitutive  equations  in  the  form 

.U>  . 1 h3rc*  E(1)  +e*  e(1)i  x(1)  . I p(1> 

*11  3 h [C11E11  +C13E33  1 » *13  " 3 h °55E13 


k(1)  ,lh3rn*  f(1)  + c*  p(1)i 
*33  3 h 1 13*11  + C33E33  1 * 


(3.36) 


where  c^,  c*3  and  c*3  are  given  in  (3.20).  Prom  (3.6)  and  (3.33)  we 


find  that  the  pertinent  plate  strains  take  the  form 
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E(1)  .1  (w*1)  +w(1))--  w(0) 

AB  2 1 A,B  B,AJ  2,AB  * 


(3.37) 


From  (3.31)  - (3.33),  we  obtain  the  usual  flexural  equations 


K(0)  +F(0)-0  K(1)  -K(0)+F(1)»O 

*B2,  B + F2  "°»  KAB,A  K2B  +PB  °’ 


(3.38) 


where  we  have  taken  proper  account  of  the  fact  that  the  thickness-shear 

stress  resultants  exist  even  though  the  associated  thickness-shear 

(0) 

plate  strains  E2A  vanish.  The  substitution  of  (3.38)2  into  (3. 38)^  yields 


+ pt1)  + P<0) 
KAB,AB  + FB,B  + F2 


(3.39) 


which  is  the  equation  of  the  flexure  of  thin  plates.  The  equation  of 
flexure  for  the  plated  crystal  plate  is  obtained  by  employing  (3.13)  in 
(3.10)  for  both  the  upper  and  lower  platings  and  then  inserting  (3.10) 
for  both  platings  in  (3.2)2,  with  the  result 


P20)  -°'  PB1>-hlKM,)A 


(3.40) 


In  this  case  of  flexure  the  displacements  w^  and  w^  , in  the  upper 


and  lower  electrode  platings,  respectively,  are  given  by 

..(0)'  (1)_  w.(0)  (0)'  _.(!)_  „(0 


(3.41) 


The  substitution  of  (3.40)  into  (3.39)  yields 

rK(1)  +h«c(0)'  K(0)#n  -0 

tKAB  +h(1CAB  KAB  ,AB  °» 


(3.42) 


Which  is  the  equation  of  static  flexure  for  the  plated  crystal  plate. 
The  substitution  of  (3.40)  into  (3.36)  yields 


^B-  0 » - C#)1,A-  4V  - 0 


(3.43) 
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which  are  consistent  with  the  integral  fonts 


where 


f "bC3*-0-  i v^3—  J “> 


(3.44) 


(3.45) 


and  A is  the  area  enclosed  by  C.  On  account  of  (3.34),  the  independent 
edge  conditions  may  not  be  obtained  directly  from  (3.44).  However,  since 
Eqs. (3.34),  (3,36),  (3.37)  and  (3.42)  - (3.44)  are  identical  in  form  with 
the  equations  of  the  elementary  theory  of  the  flexure  of  thin  plates  in 
the  absence  of  applied  loading,  the  edge  conditions  for  the  plated  plate 
are  the  same  as  in  the  elementary  theory  of  flexure  of  thin  plates.  The 
Neumann  type  uniqueness  theorem15  for  the  elementary  theory  of  the  flexure 
of  thin  plates  shows  that  the  traction  boundary  conditions  take  the  form 


where 


“aVb-5’  *£’ 


C ■ • »,vB . 


(3.46) 


(3.47) 


and  Vg  is  a unit  vector  tangent  to  c in  the  counterclockwise  direction 
and  m,  t and  v are  the  prescribed  bending  moment,  twisting  moment  and 
vertical  shearing  force,  respectively,  applied  on  an  edge.  The  alternate 
kinematic  conditions  are  on  w^0)  and  dw^/dN.  At  an  edge  of  discontin- 
uity separating  one  region  from  another,  we  require  the  continuity  of 


N»  N KL(0)  + ^ w(0) 

VaB  b * *N2  5s  \s  » 2 » 3N  2 


(3.48). 


In  this  section  we  extend  both  the  extensional  and  flexural  plate 


the  general  anisotropic  case.  However,  we  still  assume  that  the  basic 


hold 


along  with 


In  addition,  in  order  to  allow  for  the  free  thickness  strains  E, 


we  must  eliminate  the  first  order 


extensional  equation  in  (3.15)  completely  again  with  the  conditions  in 


(3.17).  From  (4.1),  (3.16)  and  (3.17),  we  have 


Since  the  temperature  change  (T-T  ) is  homogeneous,  from  (4.4) 


are  implicitly  taken 


end  (3.3)  - (3.5)  ell  first  order  piste  strains  E 


to  vanish  in  the  constitutive  equations  for  the  first  order  stress 
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resultants  even  though  the  first  order  plate  strains  E^  , Which  from  (3.6) 


are  given  by 


actually  exist  in  the  case  of  extension  by  virtue  of  the  fact  that  w^ 
and  w<2>  are  nonzero  for  extension.  If  the  first  order  plate  strains  E^^ 
are  to  be  included  in  the  constitutive  equations  for  the  the  condi- 

tions (3.16^  and  (3.17)  must  be  eliminated  and  w^  must  be  retained  in 

the  description  on  the  same  footing  as  the  wj°*  and  wJ2K  These  comments 

A A 

apply  even  in  the  case  of  the  isotropic  plate  and  pertain  to  the  descrip- 
tion developed  in  Sec. 3.  In  the  resulting  extensional  description  ob- 
tained with  the  full  inclusion  of  the  variable  at  the  very  least 

flexural  deformations  will  be  induced  in  the  general  anisotropic  case. 
These  flexural  deformations  will  not  occur  in  the  isotropic  case  or  in 
the  case  of  the  restricted  anisotropy  considered  in  Sec. 3.  The  inclusion 
of  the  in  the  constitutive  equations  for  the  and  the  attendant 

relaxation  of  conditions  imposed  and  extension  of  the  description  is 
quite  cumbersome  and  known  to  have  a negligible  influence  on  the  results 
in  the  isotropic  case  and  deemed  to  have  a sufficiently  small  influence 
on  the  results  in  the  general  anisotropic  case  that  it  may  be  omitted 
without  appreciable  error. 

Employing  (3.16)^  and  (3.17),  which  still  have  been  assumed  to  hold 
in  this  general  anisotropic  case,  in  (3.15)  we  again  obtain  (3.21)  as 
the  second  order  extensional  plate  equations.  Prom  (4.4)j_2  and  (3.3), 


we  obtain 


C’-  V 421 - Ws"  *¥  v . 
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=w°vs(Es0>  *Fi'!,s!I),V»,,-V'  (4-6> 


16 

where  we  have  introduced  the  usual  compressed  matrix  notation  for  tensor 
indices  according  to  the  scheme 

R,S  - 1,3,5  i W,V  = 2,4,6  . (4.7) 

Substituting  from  (4.6)  into  the  constitutive  equations  for  the  nonzero 
extensional  stress  resultants  in  (3.3),  we  obtain 

j>r  * t es2>]-  2,vt-  v > 


where 


Y a C _ C C-1C  ft  * V - C c-1v  . 

yrs  rs  rw  wv  vs  ’ r rw  wv  v 


(4.9) 


The  YDe  are  Voigts'  anisotropic  plate  elastic  constants  and  the  P are 

the  associated  anisotropic  plate  thermoelastic  constants.  At  this  point 

we  note  that  the  Which  exist  in  the  general  anisotropic  case, 

correspond  to  flexural  type  plate  deformations,  i.e.,  vertical  plate 

shearing  strains,  which,  however,  vanish  identically  in  ordinary  flexure. 

(2) 

The  same  sort  of  statement  holds,  of  course,  for  the  E^  . For  the  case 

of  anisotropic  extension  considered  here  the  three-dimensional  strains 

E , Which  we  need,  are  related  to  the  plate  strains  by 
KL 


Bvt  ■ -r  <WV  „ +w  ) 

KJj  < K,  Li  Li,  K. 


E 


KL 


+ X2EKL 


(4.10) 


in  Which  we  have  taken  the  first  order  plate  strains  to  vanish  as 

in  the  constitutive  equations  for  the  first  order  stress  resultants. 

The  associated  three-dimensional  displacement  fields  are  still  given  by 
(3.7).  The  plate  strains  E^^  (n«0, 2)  are  given  by 


i 


T 


E^  a i („<n>  + W<n>) 
AB  2 'a,  B B,AJ  » 


(4.11) 


and  the  remaining  plate  strains  may  be  obtained  by  solving  (4.6)  with 


the  result 


c E 
WV  VS  : 


f'VV’-V'  ^2’-- Wvs's2’-  <4-12’ 


and  we  note  that  the  relaxation  of  plate  stress  resultants  results  in 
plate  strains  and  zero  plate  rotations.  The  nth  order  plate  rotations 


take  the  form 


<n+l)  _ * (n+i ) . . 

t.  °2t  k u * 


Where 


‘Vl'I  '\k-\l1  - I 


(4.13) 


(4.14) 


As  a consequence,  we  note  for  later  use  that 

W(1)  a E(0) 

WA  E2A  » 


(4.15) 


E(0)  » i (w(0)  +w(1)>  o(0)  = 1 , M(°)  .udh  n 

The  extensional  equations  for  the  plated  crystal  plate  are  obtained 
by  employing  (3.13)  in  (3.10)  for  both  the  upper  and  lower  platings  and 
then  inserting  (3.10)  for  both  platings  in  (3.2)2  as  in  Sec. 3,  Which 
results  in  (3.25),  the  substitution  of  Which  into  (3.21)  again  yields 
(3.26)  as  the  second  order  extensional  equations  of  equilibrium  of  the, 
in  this  instance,  general  anisotropic  plate.  However,  in  this  general 
anisotropic  case  and  K™  are  given  by  (4.8),  rather  than  (3.19). 

From  (3.28)  it  is  clear  that  (3.26)  can  be  written  in  the  form 


- 0 v(2)  - 0 

*AB,A 


(4.17) 


Where,  from  (3.28),  (4.8),  (3.14),  (4.15)  and  (3.11)  for  the  upper 
plating  and  the  equivalent  of  (3.11),  but  with  primes  replaced  by  double 
primes,  for  the  lower  plating,  for  the  general  anisotropic  case  we  have 


t~H 


L 


V + 

yabcd 


C - “fan,  * T y^cd  * r + “'(I 

4 VMCD  * T '0BCD)  - “("m  '&)<*'  v - 

-2Lv'  ^(k  E(0)  +k  e (0)  - r T -r  T ) 

\h  YABCD  h YABCD/  2CBA  BA,D  2DBA  BA,C  ^20,0  fe2D  ,C  ' 

v(2)  * 2 + 3hl  v/  , 3h*  . \ (0)  2 5/ 

*AB  3 h VABCD  + h YABCD  + h YABCD/ECD  + 5 h \YABCD  + 

2L.y'  +5hlv'  V2)  -2h3(B  + V*'  +3il  V*'WT  ) - 

h YABCD  + h YABCD/ECD  3 h V AB  h VAB  h VAB/(T  0) 

h (V  VABCD  ""  "h”  YABCD)(k2CBAEBA,D+k2DBAEBA,C_  ^2CT,D~  ^2DT,C)  » <4«18) 


Where 


k * c 

ws  wv  vs  ’ 


cwvvv' 


(4.19) 


in  the  compressed  notation  and  we  note  that  T vanishes  for  the  hcno- 

> D 

geneous  temperature  states  considered  here.  The  edge  (or  boundary)  con- 
ditions at  a junction  are  still  given  by  (3.29),  (3.30)  and  the  discussion 
in  between. 

In  order  to  eliminate  extension,  Which  has  already  been  considered, 
from  the  flexural  equations,  it  is  clear  from  (3.15)  that  we  must  have 


,(0) 


(1) 


*A2 


(0) 


^2 


- 0 


(2) 


Hb 


(4.20) 


along  with 


since 


0, 


= 0. 


(4.21) 

(4.22) 


I 


In  addition,  in  order  to  allow  for  the  free  thickness  strain  , we  again 
have  (3.32).  Furthermore,  since  we  have  allowed  for  the  free  thickness 
strain  E^^,  we  must  eliminate  the  second  order  flexural  equation  in  (3.31) 
completely  again  with  the  condition  (3.33)  since 


- _ (2 ) _ 
f2  =0 


(4.23) 


Moreover,  as  usual  in  the  elementary  theory  of  flexure,  the  constitutive 
equations  for  the  zero  order  shear  stress  resultants  are  ignored  and 

the  zero  order  thickness- shear  plate  strains  are  taken  to  vanish  and 

we  again  have  (3.34).  From  (3.32)  and  (4.20)2  with  (3.3)-  (3.5)  for  m = l, 


we  obtain 


E(1)=  r-1o  F(1) 
WV  VSES  * 


(4.24) 


Where  we  have  introduced  the  conventions  shown  in  (4.7).  Substituting 
from  (4.24)  into  the  nonzero  equations  for  the  first  order  stress  result- 


ants in  (3.3),  we  obtain 


(1)  _ 2 3 (1) 

\ 3 h VRSES  » 


(4.25) 


Where  the  y are  defined  in  (4.9)  . The  remainder  of  the  equations  and 
Kb  X 

discussion  for  the  elementary  theory  of  flexure  presented  in  Sec. 3 hold 
without  change  with  the  exception  of  Eqs. (3.36)  which  are  replaced  by  (4.25) 
for  the  general  anisotropic  case  considered  in  this  section. 


5.  Large  Electrodes  on  Quartz  Plates 

A plan  view  of  the  elec trod ed  plate  is  shown  in  Fig. 2.  The  x2~ 
coordinate  axis,  Which  is  normal  to  the  major  surfaces  of  the  plate  at 
T ■ Tq,  is  arbitrarily  oriented  with  respect  to  the  crystal  axes.  Since 
the  outside  edges  of  the  plate  are  traction  free,  from  (3.29)  and  (4.17) 


we  have 


21. 


0,  *£?-0. 


(5.1) 


AB  ’ AB 

for  all  the  unelectroded  portions  of  the  plate.  Consequently,  from  (3.30) 
the  plate  edge  conditions  that  determine  the  biasing  displacement  field 
in  the  electroded  region  take  the  form 


c-°> 

V(2)  -a 
*11  “°» 

Yl0)  = 0 
*13 

*13 } = 0 ^ 

|x3l 

<X3> 

*fS’  ■ 0 - 

)r(2)  = o 

*13 

Yl0)  = 0 
*33  ° ' 

at  X3  = ±V 

|xj 

<i1.  (5.2) 

solution  satisfying 

(5.2)  and 

(4.17)  takes  the  form 

v(0> 

*AB 

(2) 

= 0 , = 0 everywhere  . 

(5.3) 

solution 

is  unique 

to  within 

static  homogeneous  plate  rotations  of 

zero  and  second  order.  Substituting  from  (4.18)  into  (5.3)  for  identical 
electrodes  on  both  surfaces,  we  obtain 

fas  * ¥ * £ yJ*<2>  - (h  * 3V  **>.  V , 


f 


Y * — 
YRS  h 


I1  . 

(5.4) 

Equations  (5.4)  constitute  six  homogeneous  linear  equations  Which  may 

readily  be  solved  for  the  six  plate  strains  eJ0^  and  Ej2\  Clearly,  if 

s s 

we  define  the  six  dimensional  vectors  and  matrices  A , B and  a a by 

or’  or  op  1 


[A^]  - [E^0),  Eg2)]  , cr-1,2  ...  6, 

<v-[(br (vir^')]-  s-1-2-6. 


i 


then  the  solution  to  (5.4)  can  be  written 

Aa-  aieVT"  v»  (5*65 

which  yields  E*;0*  and  E*2^  as  known  expressions  linear  in  (T-T  ).  Note 

bb  O 

that  the  solution  is  independent  of  f and  l and  further,  that  since 
(2) 

and  constitute  two  planar  tensors  that  vanish  everywhere,  the  plate 


edge  conditions 


naC  * °>  naC  - °- 


at  the  edge  of  an  electrode  are  satisfied  for  the  edge  having  a variable  N. 

A 

and,  hence,  being  a curve.  Consequently,  the  solution  is  valid  for  elec- 
trodes of  arbitrary  shape  provided  the  electrodes  are  large  compared  to 
the  thickness  of  the  quartz  plate*8.  When  E^°^  and  E^2^  have  been  de- 
termined  from  (5.6),  E^°^  and  E^2^  are  readily  determined  from  (4.12). 
lhen  the  three-dimensional  strain  can  be  evaluated  from  (4.10),  Which 
holds  for  this  purely  extensional  case.  There  is  no  flexure  because  the 

temperature  change  (T-T  ) is  homogeneous  and  the  electrodes  are  identical. 

o 

It  is  well  known  that  in  static  linear  elasticity  the  solution  to  a 

boundary  value  problem  is  unique  only  to  within  a static  homogeneous 

19 

(global)  infinitesimal  rigid  rotation  . In  addition,  the  change  in  fre- 
quency due  to  a homogeneous  infinitesimal  rigid  rotation  is  shown  to 
vanish  in  the  Appendix.  Consequently,  without  any  loss  in  generality, 
we  may  select  the  homogeneous  rigid  rotation  to  take  any  value  that  is 
convenient  and  in  particular  to  vanish.  Accordingly,  we  take 
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which  with  (4.10)  yields 

W =3  E 

L,K  KL 


3(0) 

KL 


* 4-S’ 


(5.8) 


(5.9) 


which  with  (5.6)  provides  the  biasing  displacement  gradients  w as  a 

N 

known  linear  function  of  {■£  “ TQ) • Thus,  we  may  now  substitute  from  (5.6) 

into  (4.12)  and  then  from  (5.6)  and  (4.12)  into  (5.9),  Which  may  then  be 

substituted  into  (2.11)  with  (2.12)  to  obtain  c as  a known  linear 

Lynof  — — . 

function  of  (T  - T ) . 

o 


6.  Thickness  Modes  in  Piezoelectric  Plates 

Thickness  modes  in  piezoelectric  plates  are  solutions  depending  on 
the  thickness  coordinate  only,  which  satisfy  the  linear  piezoelectric 
equations  and  the  boundary  conditions  on  the  major  surfaces  of  the  plate, 
but  do  not  satisfy  any  conditions  on  the  minor  surfaces.  The  thickness 
solutions  are  of  practical  importance  because  plates  with  lateral  dimen- 
sions very  large  compared  to  the  thickness  are  employed  in  resonators  and 
filters,  and,  consequently,  the  actual  mode  in  the  bounded  plate  is  not 
that  different  from  the  thickness  mode.  In  any  event  in  this  first  treat- 
ment of  temperature  induced  frequency  changes  in  electroded  quartz  plates 
only  the  thickness  modes  will  be  considered. 

Since  is  the  thickness  coordinate,  we  substitute  from  (2.3)  into 
(2.4)  and  retain  Xj-dependence  only  to  obtain 


22y«2  «,  22 


W, 


22  “ PV 


-0 


V 


e22y\f,  22  “ *22^,22 


(6.1) 
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Which  are  the  differential  equations  that  must  be  satisfied  by  the  thick- 
ness solution.  For  the  thickness  eigensolutions  of  interest  here,  from 
(2.14)  for  the  linear  case  only,  with  (2.3)^  and  retaining  X2~dependence 
only,  we  obtain  the  boundary  conditions 


22yor2Ua,  2 +e22y^,  2 * ^ 2h  p atx2“±h. 


(6.2) 


since  the  electrodes  are  shorted.  The  solution  for  thickness-modes  in 

arbitrarily  anisotropic  piezoelectric  plates  with  shorted  electrodes  on 
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the  major  surfaces  may  be  written  in  the  form 

3 

V n(n)R(n)  v v 

l B PY  81"  W » 


\ 


n=l 
3 


?-[t?  i ’,n>>r 

n n»l 


(6.3) 


which  satisfies  (6.1)  provided 

- (n) 


(C2yof2  " C 


, (n) 


6 )0(n) 
yof  a 


. -(n)  2 2 

0 , c = paj  /T^  , 


(6.4) 


and  the  0^“'  are  the  normalized  eigenvectors  of  the  linear  homogeneous 


algebraic  system  in  (6.4)  for  the  eigenvalues  c 
stiffened  elastic  constants,  which  are  given  by 


- (n) 


of  the  piezoelectrically 


C2yor2  " 22ya2  + ®22ye22c/e22  * 

In  order  that  (6.2)2  be  satisfied  we  must  have 


and  (6.2)^  are  satisfied  if  the  are  given  by 

B<n».-  6^’.32t(V5'n,\l»os  \h-R\h  sin  T^hl  , 


(6.5) 


(6.6) 


(6.7) 


where  we  have  employed  (6.4)  and 
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pWpfoO  _ 6^,  R “ 2p/h//ph, 


(6.8) 


in  obtaining  (6.7)  from  (6,2)^.  Substituting  from  (6.7)  into  (6.6),  we 
obtain  3 

ifi-  Y (k(n))2/T^h(cot  T^h-RT^h)^]  = 0,  (6.9) 


where 


n=l 


q (n)  _ B(n) 

(k(n)  2 _ py  e22ypor  e22g 

5(n)e 
C ®22 


(6.10) 


The  condition  for  a nontrivial  solution  of  the  scalar  equation  (6.9)  is 
3 

Y <* (n>  > 2/\h  (cot  T\k  - R71nh ) = 1 , (6.11) 

n«l 


the  roots  of  which  determine  the  resonant  frequencies  of  thickness  vibra- 
tion of  piezoelectric  plates  driven  by  the  application  of  a voltage  across 
the  surface  electrodes.  From  this  unperturbed  thickness  eigensolution 
we  can  determine  the  normalized  eigensolution  we  need  for  the  perturba- 
tion formulation  in  Sec. 2 simply  by  writing 


<3y  - vya , f-$/K, 


Where,  from  (2.8),  including  the  electrode  platings, 
we  find 


»j-p  y 


3 

I 


3 


I 


« (m)o  (m)_ (n)D  (n)T  Sin(T1m"  Vh 

Y ^ L — ’ 

+ 2Rh  sin  l^h  sin  T^hJ  . 


(6.12) 

(6.3)1  and  (6.8>2 
sin  (1^  + 11  )h 

' vC 

(6.13) 


26. 


7.  Temperature  Dependence  of  Resonant  Frequency 

Ttoe  change  in  the  frequency  of  thickness  vibrations  with  temperature 
of  any  electroded  quartz  plate  may  now  be  determined  from  (2.1),  which  we 
rewrite  here  for  any  one  mode  in  the  form 


am“V2V  ^ 


(7.1) 
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where  for  the  case  of  thickness  vibrations  we  have 

h 


»H--  I ^<,2  aV 


(7.2) 


-h 


From  (2.5)  and  (2.10)  for  the  thickness-mode  being  perturbed  here,  we  have 

M 


*2y  “ (c2y2er  + A22y2cr)ga,  2 » 


(7.3) 


where  ^2y2a  (2. 11), with  (2. 12), is  known  as  a linear  expression 

in  (T-T  ) from  (5.9), with  (5.5),  (5.6)  and  (4.12)  and 


(7.4) 


The  d^2y2o^dT  are  obtained  from  the  first  temperature  derivatives  of  the 

6 ~ 


fundamental  elastic  constants  of  quartz  dc^^^/dT  referred  to  the  principal 
axes  by  the  tensor  transformation  relation 


d ^ d ~ 

dT  §2y2er  " a2D*YEa2FaoG  dT  f DEFG  * 


(7.5) 


where  the  a are  the  matrix  of  direction  cosines  for  the  transformation 

yt, 


from  the  principal  axes  to  the  coordinate  system  containing  the  axis  normal 
to  the  plane  of  the  plate.  When  the  conventional  IEEE  notation22  for 
doubly-rotated  plates  is  written  in  the  form  (Y,X,w, /)<p,  0,  where  ♦ ■ 0,  the 
rotation  angles  <p  and  0 are  the  first  two  Euler  angles,  and  the  a^  are 
given  by 


s 

- 
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cos  q>  sin  cp  0 

- cos  0 sin  cp  cos  6 cos  cp  sin  6 

sin  6 sin  cp  - sin  9 cos  cp  cos  6 


_L  > 


(7.6) 


since  f - 0.  Clearly,  the  transformation  relations  for  the  second  and 


third  order  elastic,  piezoelectric  and  dielectric  constants,  and  coeffi- 
cients of  linear  expansion  may  be  written  in  the  respective  forms 


jKLMN  ° aKD3LEaMFaNG2DEFG  * 

<V 

SkLMNAB  = aKDaLE  MF  NGaAHaBI§DEFGHI  * 

e & a a a e 
KLM  KD  LE  MF  DEF 

e = a a e . ot  — a a a 
KL  KM  LN  MN  ’ KL  KM  LN  MN * 


(7.7) 


Where  the  tensor  quantities  with  the  upper  cycle  are  referred  to  the 
principal  axes  of  the  crystal. 

Calculations  of  the  resonant  frequency,  called  ou  here,  for  the 
particular  unperturbed  thickness-mode  of  interest  proceed  by  numerically 
finding  the  value  of  <uu  satisfying  (6.4)  and  (6.11).  After  the  cu  of 
interest  has  been  determined,  the  full  Mth  eigensolution  is  obtained  by 
substituting  the  and  T)^  from  (6.4)  and  the  from  (6.7)  into  (6.3). 

The  normalization  integral  in  (6.13)  is  then  evaluated  and  the  normalized 
eigensolution  obtained  from  (6.12).  Then  the  perturbation  integral  is 
evaluated  by  employing  (7.3)  and  (6.12)^  in  (7.2)  and  performing  the 
integrations.  On  account  of  (5.9),  with  (5.5),  (5.6),  (4.12)  and  (6.3)^  two 
distinct  types  of  integrals  arise  in  (7.2),  the  evaluation  of  Which  yields 


n 

P M M 

J 9Qf,29Y, 


...  IV  V r(p)p(p)-(q)D(q)T1  r®in  (Wh 

2^2-7  l 1 B B PY  W <W  ' 

nm  p=l  q=l  , P 4 

" si»  <VVhi 

* 01  -V)  J ’ 

p q 


3 3 


P —2 _ M M ...  1 v V „(p)D(p)_(q)0(q)^  *.  17  h 

J*aW“a-7  il1  B 8V  WVv^ 


p=i  q=i 


<VV 


r)  sin(T>P+Vh  + 


% + \ 


'<VVh  + ^ir^T ' 


p q 


(YV 


•)  sintTlp-Tl^h  + 2k--  2 cos (^ - T]g)h]  . (7.8) 


Calculations  have  been  performed  using  the  known  values  of  the  second 


order  elastic,  piezoelectric  and  dielectric  constants  of  quartz  , the 
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third  order  elastic  and  thermoelastic  constants  of  quartz  and  the 


recently  obtained  * temperature  derivatives  of  the  fundamental  elastic 


constants  of  quartz.  Specific  calculations  have  been  made  for  doubly- 


rotated  quartz  plates  1.7  mm  thick,  with  4000  A thick  gold  electrodes  on 


the  major  surfaces,  vibrating  in  the  fundamental  thickness-mode.  Typical 


results  of  the  calculations  are  shown  in  Figs.3-8.  Figure  3 shows  the 


actual  change  in  frequency  for  an  electroded  quartz  plate  for  a fixed 


value  of  9 - - 49.2167°,  Which  corresponds  to  the  BT-cut  for  cp  - 0 , as  a 


function  of  <p.  The  lower  and  upper  curves  are  for  the  B-  and  C-modes, 


respectively,  Which  are  defined  so  that  for  each  cut  and  mode  number  n 


the  sequence  f>  f_>  f is  followed.  The  curves  in  the  figures  repeat 

ABC 


because  of  the  synmetry  of  quartz.  Figure  4 shows  the  change  in  frequency 


due  to  the  presence  of  the  electrodes  for  the  same  orientations  and  modes 


of  the  electroded  plate  shown  in  Fig. 3.  Note  that  the  change  in  frequency 


due  to  the  electrodes  is  about  two  orders  of  magnitude  smaller  than  the 


i 


fa 


r 


i 
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actual  change  shown  in  Fig. 3.  Since  a portion  of  the  aging  rate  is  a 
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result  of  the  relaxation  of  residual  stress  in  the  electrodes  , the 
change  in  frequency  shown  in  Fig. 4 is  the  portion  of  the  actual  change 
in  frequency  that  contributes  to  the  aging  rate.  Figure  5 shows  the  actual 
change  in  frequency  for  an  electroded  quartz  plate  for  a fixed  value  of 
9 - 35.25°,  Which  corresponds  to  the  AT-cut  for  9*0°,  as  a function  of  9. 
Note  that  the  C-mode,  which  is  the  piezoelectrically  active  mode  for  the 
AT-cut,  i.e.,  for  9=0°,  is  relatively  flat  and  has  no  zero  crossings, 
while  the  B-mode  has  wide  excursions  and  zero  crossings.  Hie  zero  cross- 
ings for  the  B-mode  correspond  closely  to  points  on  the  locus  of  zero 
temperature  coefficients  shown  in  Fig. 2 of  Ref. 26.  Figure  6 shows  the 
change  in  frequency  due  to  the  presence  of  the  electrodes  for  the  cases 
for  which  the  actual  changes  in  frequency  are  shown  in  Fig. 5.  The  dotted 
line  in  Fig. 7 represents  the  actual  change  in  frequency  near  a zero 
crossing  to  an  enlarged  scale  for  a fixed  value  of  9=5°  for  the  B-mode, 
and  the  solid  line  represents  the  change  in  frequency  due  to  the  presence 
of  the  electrodes.  Note  that  the  two  lines  have  different  scales,  one 
being  two  orders  of  magnitude  smaller  than  the  other.  The  difference  in  6 
for  the  zero  crossings  is  related  to  what  is  called  the  apparent  shift  in 
angle  of  the  zero  temperature  cut  with  electrode  thickness.  Hie  vertical 
distance  from  the  intersection  of  the  dotted  line  with  the  horizontal  axis, 
which  is  the  actual  zero  temperature  cut  of  the  1.7  mm  thick  quartz  plate 
with  4000  A thick  electrodes,  to  the  solid  line  measured  on  the  inner  scale 
represents  the  change  in  frequency  that  can  contribute  to  aging.  Thus, it 
is  desirable  to  find  actual  zero  temperature  cuts  which  have  that  ordi- 


nate as 


ill  as  possible  for  a given  electrode  thickness.  Hie  dotted 
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line  in  Fig. 8 represents  the  actual  change  in  frequency  near  a zero  cross- 
ing to  an  enlarged  scale  for  a fixed  value  of  cp  - 30°  for  the  C-mode,  and 
the  solid  line  represents  the  change  in  frequency  due  to  the  presence  of 
the  electrodes.  The  zero  crossings  in  both  Figs. 7 and  8 correspond  closely 
to  points  on  the  locus  of  zero  temperature  coefficients  shown  in  Fig. 2 of 
Ref. 26,  as  they  should. 

Acknowledgements 

We  wish  to  thank  D.  Stevens  of  Rensselaer  Polytechnic  Institute  for 
help  with  the  calculations. 

This  work  was  supported  in  part  by  the  Army  Research  Office  under 
Grant  HO.  DAAG29-76-G-0173,  the  Office  of  Haval  Research  under  Contract 
No.  N00014-76-C-0368  and  the  National  Science  Foundation  under  Grant 
No.  ENG  72-04223. 


31 


REFERENCES 


1.  H.F.  Tiersten,  "Perturbation  theory  for  Linear  Electroelastic  Equations 
for  Small  Fields  Superposed  on  a Bias,"  j.  Ac oust.  Soc.  Am.,  64, 

832  (1978). 

2.  B.K.  Sinha  and  H.F.  Tiersten,  "On  the  Influence  of  a Flexural  Biasing 
State  on  the  Velocity  of  Piezoelectric  Surface  Waves,"  Wave  Motion,  1, 
37  (1979). 

3.  R. D.  Mindlin,  "An  Introduction  to  the  Mathematical  Theory  of  the 
Vibration  of  Elastic  Plates,"  U. S.  Army  Signal  Corps  Eng.  Lab., 

Fort  Monmouth,  New  Jersey  (1955).  Signal  Corps  Contract  DA-36-039SC- 
56772. 

4.  R.D.  Mindlin,  "High  Frequency  vibrations  of  Crystal  Plates,"  Quart. 
Appl.  Math.,  19,  51  (1961). 


5. 

6. 


7. 

8. 


9. 


10. 

11. 


12. 

13. 

14. 


H.F.  Tiersten,  Linear  Piezoelectric  Plate  Vibrations  (Plenum,  New 
York,  1969),  Chap. 13. 

B.K.  Sinha  and  H.F.  Tiersten,  "First  Temperature  Derivatives  of  the 
Fundamental  Elastic  Constants  of  Quartz,"  to  be  published  in  the 
J.  Appl.  Phys. 

Ref. 5,  Chap. 14,  Sec. 4. 

J.  Tasi  and  G.  Herrmann,  "Thermoelastic  Dissipation  in  High  Frequency 
Vibrations  of  Crystal  Plates,"  J.  Acoust.  Soc.  Am.,  3£,  100  (1964). 

R.D.  Mindlin,  "High  Frequency  Vibrations  of  Plated  Crystal  Plates," 
in  Progress  in  Applied  Mechanics  (Macmillan,  New  York,  1963), 
pp. 73-84. 

This  procedure,  which  follows  the  work  of  Mindlin,  results  in  a 
valuable  reduction  and  simplification  in  the  equations. 

Since  does  not  vanish,  the  conditions  in  Eq.  (3.17)  are  implied 

by  (3.15)2  and  (3.16)^.  Certain  consequences  of  (3.17)  are  discussed 
in  detail  in  Sec. 4. 

A.E.H.  Love,  A Treatise  on  the  Mathematical  Theory  of  Elasticity. 

4th  ed.  (Cambridge  University  Press,  Cambridge^  1927 ) also  (Dover, 

New  York,  1944)  Sec. 118. 

Ref. 5,  Chap. 13,  Sec. 3. 

(2) 

Since  w2  does  not  vanish,  the  conditions  in  (3.33)  are  implied  by 
(3.31)3  and  (3.32).  The  second  order  plate  equations  are  not  im- 
portant in  the  elementary  theory  of  the  flexure  of  thin  plates. 


■fek 


32 


15.  Ref. 3,  Sec. 6. 04. 

16.  Ref. 5,  Chap. 7,  Sec.l. 

17.  The  nth  order  plate  rotations  are  obtained  by  substituting  from  the 
thickness  expansion  for  the  mechanical  displacement  field  into  the 
expression  for  the  three-dimensional  local  small  rotation  and  pro- 
ceeding exactly  as  in  the  determination  of  the  nth  order  plate 
strains  in  (3.6). 

18.  Although  the  approximate  plate  equations  yield  accurate  results  for 
points  not  too  close  to  the  edges  of  the  electrodes,  the  results 
are  not  at  all  meaningful  in  the  near  vicinity  of  the  electrode 
edges.  This  is  a result  of  the  fact  that  the  three-dimensional 
elastic  strain  field  has  a singularity  along  the  edge  of  the  elec- 
trode. However,  even  in  the  least  ductile  of  real  materials,  the 
inelastic  behavior  prevents  the  singularity  from  occurring  in  the 
actual  stress  and  strain  distributions,  which  remain  finite.  The 
stress  and  strain  distributions  in  the  vicinity  of  the  electrode 
edges  obtained  from  the  approximate  plate  equations  are  quite  dif- 
ferent from  the  actual  ones  that  occur  even  when  the  inelastic 
behavior  is  taken  into  account.  Nevertheless,  this  difference  is 
not  important  for  our  purposes  because  in  the  case  of  large  area 
electrodes,  in  Which  we  are  interested,  the  strain  distribution 
obtained  from  the  approximate  plate  equations  is  accurate  over 
almost  the  entire  electroded  region,  and  the  small  area  in  the 
vicinity  of  the  electrode  edges  does  not  have  a significant  influence. 

19.  Ref. 12,  Secs. 18  and  118. 

20.  H.F.  Tiersten,  "Thickness  Vibrations  of  Piezoelectric  Plates," 

J.  Acoust.  Soc.  Am.,  35,  53  (1963). 

21.  For  a wide  plate  with  wide  electrodes  the  influence  of  energy 
trapping  can  be  ignored. 

22.  "Standards  on  Piezoelectric  Crystals,  1949,"  Proc.  IRE,  37^  1378 
(1949). 

23.  R.  Bechmann,  "Elastic  and  Piezoelectric  Constants  of  Alpha-Quartz,” 

Phys.  Rev.,  1060  (1958). 

24.  R.N.  Thurston,  H. J.  McSkimin  and  P.  Andreatch,  Jr.,  "Tfoird  Order 
Elastic  Constants  of  Oiartz,"  J.  Appl.  Phys.,  37,  267  (1966). 

25.  F.  Kohlrausch,  Lehrbuch  der  prakt.  Fhysik,  16.  aufl.  5.  158  (1930). 

26.  R.  Bechmann,  A.D.  Ballato  and  T. J.  Lukaszek,  "Higher  Order  Temperature 
Coefficients  of  the  Elastic  Stiffnesses  and  Compliances  of  Alpha 
Quartz,"  Proc.  IRE,  50,  1812  (1962). 


27.  E.  P.  Eer  Nisse,  "Simultaneous  Thin-Film  Stress  and  Mass-Change 
Measurements  Using  Quartz  Resonators,"  J.  Appl.  Ifcys. , 43,  1330 
(1972). 

28.  Ref. 5,  Chap. 8,  Sec. 7.  In  this  instance  the  operations  are  performed 
for  more  than  one  region. 

29.  The  sign  of  the  second  term  in  Eq. (2.12)2  of  Ref . 1 should  be  changed 
from  a minus  to  a plus. 


mtmm 


S 


Figure  1 
Figure  2 
Figure  3 

I 

Figure  4 

Figure  5 

i 

j 

Figure  6 

Figure  7 


Figure  8 


FIGURE  CAPTIONS 


Schematic  diagram  of  the  plated  crystal  plate 

Plan  view  of  rectangular  electrode  on  quartz  plate 

Actual  relative  change  in  the  fundamental  resonant 
frequency  per  °K  for  1.7  mm  thick  doubly-rotated  quartz 
plate  with  4000  A thick  gold  electrodes  as  a function 
of  cp  with  0--  49°  13' . 

Relative  change  in  the  fundamental  resonant  frequency 
per  °K  due  to  the  electrodes  for  the  cases  treated  in 
Fig. 3. 

Actual  relative  change  in  the  fundamental  resonant 
frequency  per  °K#for  a 1.7  mm  thick  doubly-rotated  quartz 
plate  with  4000  A thick  gold  electrodes  as  a function 
of  9 with  0-35°  15' . 

Relative  change  in  the  fundamental  resonant  frequency 
per  °K  due  to  the  electrodes  for  the  cases  treated  in 
Fig.  5 

Relative  change  in  the  fundamental  resonant  frequency 
per  ° K for  a 1.7  mm  thick  doubly-rotated  quartz  plate 
with  4000  A thick  gold  electrodes  near  a zero  crossing  as 
a function  of  0 with  9 ■ 5° . The  dotted  and  solid  lines 
indicate  the  actual  and  electrode  induced  changes  in 
frequency,  respectively. 

Relative  change  in  the  fundamental  resonant  frequency 
per  °K  for^a  1.7  mm  thick  doubly-rotated  quartz  plate 
with  4000  Athick  gold  electrodes  near  a zero  crossing 
as  a function  of  0 with  9«30°.  The  dotted  and  solid 
lines  indicate  the  actual  and  electrode  induced  changes 
in  frequency,  respectively. 
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APPENDIX 


We  now  show  from  the  perturbation  integral  for  A^  in  (2.1)  and  (2.2) 
that  A^  vanishes  for  an  arbitrary  pure  homogeneous  infinitesimal  rigid 
rotation  Q^,  which  is  given  by 

<Al> 

Since  under  these  circumstances  the  strain  vanishes  as  does  the 

KL 

stress  T^,  from  (2.12)2  and  (Al)  we  have 


^KL  * °Kl"  ” ^LK’ 


and  from  (2.11)  we  obtain 


c ■ c n1  + c 

LYMU  2 LY KM  KV  2LKMV  Ky 


Since  the  biasing  displacement  field  resulting  from  the  nonlinear 
behavior  of  the  electrode  plating  has  not  been  included  in  the  description, 
i.e.,  in  Eq. (2.14)^,  the  form  of  the  perturbation  integral  in  (2.2)  does 
not  actually  result  in  a zero  A^  for  nonzero  oj^.  However,  if  the  bias 
due  to  the  nonlinear  behavior  of  the  electrode  plating  is  properly  in- 
cluded in  the  description,  the  amended  form  of  the  perturbation  integral 
does  result  in  a zero  ^ for  nonzero  0^.  The  unbiased  plating  equations, 
Which  enable  the  entire  effect  of  the  plating  to  be  treated  as  a homogene- 
ous boundary  condition  at  the  surface  of  the  piezoelectric  plate  and 
thereby  result  in  a major  simplification  in  the  analysis,  are  based  on  a 
number  of  simplifying  thin  plate  assumptions7  which  result  in  the  occur- 
rence of  Voigts'  anisotropic  linear  plate  constants.  Since  the  use  of 
the  approximate  thin  plate  assumptions  in  the  case  of  the  biased  plating 
will  result  in  the  occurrence  of  effective  anisotropic  plate  constants 
other  than  the  Voigt  constants,  it  does  not  appear  to  be  particularly 


A2 


purposeful  to  obtain  the  biased  approximate  thin  plating  equations.  Con- 
sequently, in  this  appendix  we  take  the  alternative  course  of  extending 
the  perturbation  integral  for  to  the  case  Where  the  electrode  platings 
are  included  as  additional  three-dimensional  regions  attached  to  the 
piezoelectric  plate.  We  then  show  that  the  extended  perturbation  integral 
results  in  a zero  for  nonzero  Q^.  However,  in  order  to  keep  the  demon- 
stration simple  and  clear  and  not  introduce  additional  complications  we 
first  show  that,  for  no  electrode  plating  on  a purely  elastic  plate,  the 
perturbation  integral  presently  in  (2.18)  vanishes  for  an  arbitrary  , 
and  we  further  note  that  since  the  actual  electrode  plating  is  very  thin, 
the  perturbation  integral  in  (2.18)  very  nearly  vanishes  for  thin  electrode 
platings  on  the  piezoelectric  plate. 

Substituting  from  (A3)  into  (2.19)  for  zero  temperature  change,  we 

obtain 
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which  by  virtue  of  the  symmetry  of  the  ^LJKM  and  the  fact  that  0^ 
homogeneous,  may  be  written 


is 
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Substituting  from  (2.9),  for  the  purely  elastic  case  (e  * 0),  we  obtain 
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Which  with  the  aid  of  the  divergence  theorem  and  the  normalized  eigen- 
solution  form  of  (2.4)^  yields 
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Since  for  traction-free  boundary  conditions  in  the  linear  eigenvibration 
problem  we  have 
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(A8) 


Eqs. (2.1),  (A2)2  and  (A7)  yield 
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In  order  to  obtain  the  appropriate  form  of  the  perturbation  integral 
for  the  electroded  piezoelectric  plate,  we  return  to  Eqs. (3.17)  and  (3.18) 
of  Sef.l,  Which  have  the  form 
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where  (A10)  are  the  equations  satisfied  by  the  ;tth  eigensolution  and 
(All)  are  the  equations  satisfied  by  the  nearby  perturbed  solution  at 
frequency  <u.  To  (A10)  and  (All)  we  must  adjoin  the  equations  for  the 
nonpiezoelectric  electrodes.  Which  may  be  written  in  the  form 


bY,b 


(A12) 

5^.0, 

(A13) 

Where  m represents  the  mth  electrode,  say  m •»  1 for  the  top  electrode  and 
m » 2 for  the  bottom  electrode,  and  since  we  are  concerned  with  perturba- 
tions in  Which  the  electrodes  remain  shorted,  we  have  (A12)2  and  (A13)2> 
From  the  first  of  the  conditions  in  (A10)  - (A13),  we  form 
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and  at  this  point  for  our  purposes,  from  (2. 3 and  (2.5)^,  we  note  that 
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where  from  (2.11)  for  zero  biasing  stress  we  have 
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Performing  the  usual  operations  , employing  (2.3),  (A10)2,  (All)2  and 
the  divergence  theorem  in  the  usual  manner,  we  obtain 
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Since  the  perturbed  solution  is  nearby  the  unperturbed  solution,  we  have 

|A|  «®i,  IHyl  (A18) 

Substituting  from  (A18)  into  (A17),  neglecting  products  of  small  quantities, 
introducing  the  normalization  integral  and  employing  (2.3)  for  the  lith  mode, 
(2.7)  and  (2.9),  we  obtain 
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In  (A20)  the  integrals  over  S are  of  two  types,  one  in  which  the 
surface  abuts  free-space  and  the  other  in  Which  the  surface  abuts  the 
piezoelectric  body.  Similarly,  since  the  piezoelectric  body  is  only 
partially  electroded,  the  integral  over  is  of  two  types  also,  one  in 
which  the  surface  abuts  free-space  and  the  other  in  which  the  surface 
abuts  an  electrode.  At  the  interfaces  between  the  electrodes  and  the 
insulator,  the  surface  element  in  (A20)  occurs  twice,  once  from  each  side 
with  outwardly  directed  unit  normals  N , which  then  are  oppositely  di- 
rected. If  we  agree  to  count  each  such  surface  only  once  and  adopt  the 
convention  that  the  normal  N.  is  positive  when  directed  out  of  the  insula- 

Jj 

tor  into  the  electrode  at  a contiguous  surface,  (A20)  can  be  written 
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in  Which  * denotes  the  portion  of  the  surface  of  the  insulator  abutting 
free  space,  denotes  the  portion  of  the  surface  of  the  insulator  abutting 
an  electrode  and  denotes  the  portion  of  the  surface  of  the  electrodes 
abutting  free-space,  and  Where  we  have  employed  the  conditions 

^ m <§,  g"  - gy  , <A24> 

along  the  surfaces  on  which  the  electrodes  are  attached  to  the  insulator. 
Since  the  nth  piezoelectric  eigensolution  satisfies  the  boundary  conditions 
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from  (A23)  H can  be  written  in  the  form 
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(A26) 
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In  the  perturbation  integral  in  (A26 ) the  quantities  without  the  n 
are  perturbation  terms  Which  are  to  be  determined  from  the  nth  eigensolu- 
tion  due  to  the  presence  of  the  bias.  The  boundary  conditions  that  exist 
in  the  presence  of  the  bias  are 
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Substituting  from  (A27)  into  (A26 ) , rewriting  the  surface  integrals  so 
that  normals  to  all  surfaces  are  positive  as  in  Eq. (A20)  and  employing  the 
divergence  theorem,  we  obtain 
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which  is  the  form  of  the  perturbation  integral  of  particular  interest  to 
us.  In  (A28)  the  variables  with  the  superscript  n take  the  values  given 
by  the  nth  orthonormal  eigensolution  g^,  f^,  g^  in  the  presence  of  the 
bias,  and  consequently,  from  (2.5)  for  zero  temperature  change,  (A15) 2 
and  (A21)  we  have 
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Where  from  (2.12)  of  Ref.l,  including  the  terms  depending  on  the  biasing 

deformation  only.  Which  are  the  only  ones  needed  for  our  purposes  here, 
29 


we  have 
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where  eQ  ant^  denote  the  first  order  electroelastic  constants, 

the  electric  permittivity  of  free-space  and  the  Jacobian  of  the  biasing 
deformation,  which  is  approximately  unity.  The  substitution  of  (A29) 
into  (A28)  yields 
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Equation  (A31)  is  the  form  of  the  perturbation  integral  we  need  to  show 
that  Att^  vanishes  for  a nonzero  pure  homogeneous  infinitesimal  rigid 
rotation  0^  of  the  electroded  piezoelectric  body. 

For  a homogeneous  infinitesimal  rigid  rotation  of  the  electroded 
piezoelectric  body,  from  (A2),  (A16)  and  (A30),  we  have 
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along  with  (A3),  since  E vanishes  and  the  electrodes  are  attached  to 

Kli 

the  piezoelectric  body.  Substituting  from  (A3)  and  (A32)  into  (A31),  we 
obtain 
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Which  by  virtue  of  the  synmetry  of  the  c , e„„  and  c_  and  the  fact  that 
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is  homogeneous,  may  be  written 
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Substituting  from  (2. 9)  ^ and  (A15),  with  (A21)  and  omitting  the  super- 


script t,  we  obtain 
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which  with  the  aid  of  the  divergence  theorem,  the  normalized  forms  of  (AlO)^ 
and  (A12>1  and  the  conventions  introduced  in  (A23)  yields 
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where  we  have  employed  (A24)  in  obtaining  (A36).  Now,  from  (A19)  and  (A36), 
with  the  conditions  in  (A2)2  and  (A25),  we  have 

^ = 0 . (A37) 


